ASYMPTOTICS OF CLASSICAL SPIN NETWORKS 



STAVROS GAROUFALIDIS AND ROLAND VAN DER VEEN 

Abstract. A spin network consists of a cubic ribbon graph (i.e., an abstract trivalent graph with a cychc 
ordering of the vertices around each edge) and an admissible coloring of its edges by natural numbers. The 
evaluation of a spin network, with our normalization, is an integer number. The main results of our paper 
are: (a) an existence theorem for the asymptotics of evaluations of arbitrary spin networks (using the theory 
of G-functions) , (b) a rationality property of the generating series of all evaluations with a fixed underlying 
graph (using the combinatorics of the chromatic evaluation of a spin network), (c) a complete study of the 
asymptotics of the so-called 6j-symbols in all cases: Euclidean, Plane or Minkowskian (using the theory of 
Borel transform), (d) an assignement of two numerical invariants (a spectral radius and a number field) to 
spin networks and cubic graphs and (e) an explicit illustration of our results for some Euclidean, Plane and 
Minkowskian 6j-symbols (using the constructive Wilf-Zeilberger method). 
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1. Introduction 

1.1. Spin networks and their evaluations. The paper is concerned with 

(a) an existence theorem for the asymptotics of evaluations of arbitrary spin networks (using the theory 
of G- functions), 

(b) a rationality property of the generating series of all evaluations with a fixed underlying graph (using 
the combinatorics of the chromatic evaluation of a spin network) , 

(c) a complete study of the asymptotics of the so-called 6j-symbols in all cases: Euclidean, Plane or 
Minkowskian (using the theory of Borel transform) , 

(d) the assignement of two numerical invariants to spin networks and cubic graphs, the spectral radius 
and a number field (see Definition 12. 3p , 

(e) an explicit illustration of our results for the regular and some flat 6j-symbols (using the constructive 
WZ method). 

A (classical) spin network (r,7) consists of a cubic ribbon graph T (i.e., an abstract trivalent graph with 
a cyclic ordering of the edges at each vertex) and a coloring 7 of its set of edges by natural numbers. 

According to Penrose, spin networks correspond to a diagrammatic description of tensors of representations 
of SU(2), and their evaluation is a contraction of the above tensors. Spin networks originated in Racah's 
work in atomic spectroscopy in the late forties \Ka\ and in the work of Wigner in the fifties [Wi] and their 
evaluations (exact or asymptotic) is a useful and interesting topic studied by Ponzano-Regge, Biedenharn- 
Louck and many others; see [BLlllBL2llPRirV"MKj . In the past three decades, spin networks have been used 
in relation to classical and quantum gravity and angular momentum in 3-dimensions; see |EPRi [PeTl IPe2i [RS] . 
On the mathematical side, q-deformations of spin networks (so called quantum spin networks) appeared in 
the eighties in the work of Kirillov-Reshetikhin [KR] . Quantum spin networks are knotted framed trivalent 
graphs embedded in 3-space with a cyclic ordering of the edges near every vertex, and their evaluations are 
rational functions of a variable q. The quantum 3j and 6j-symbols are the building blocks for topological 
invariants of closed 3-manifolds in the work of Turaev-Viro |TV1 ITu] . Quantum spin networks are closely 
related to a famous invariant of knotted 3-dimensional objects, the celebrated Jones polynomial, [j]. A 
thorough discussion of quantum spin networks and their relation to the Jones polynomial and the Kauffman 
bracket is given in [KLj and [CFS| . 

Aside from the above applications, spin networks and their evaluations have an intrinsic interest, they are 
formulated in terms of elementary combinatorial data and offer many challenges even for relatively simple 
networks such as the cube. Some examples of spin networks that will be discussed in the paper are shown 
in Figure [TJ 




Figure 1. From left to right: The theta or 3j-symbol, the tetrahedron or 6j-symbol, the Cube, the 5-sided 
drum and the complete bipartite graph ^"3,3 or 9j-symbol. The cyclic order of the edges around each vertex 
is counterclockwise. 



Spin networks with underlying graph or a tetrahedron are usually called 3 j -symbols and 6 j -symbols 
respectively and they form the building blocks for the evaluation of any spin network; see for example Section 
l6.1l or [BLl|,[BL2llPRj . We now discuss the evaluation of a spin network, introduced by Penrose in [Pel|,IPe2] . 

Definition 1.1. (a) We say a spin network is admissible when the sum of the three labels a,b,c around 
every vertex is even and a, b, c satisfy the triangle inequalities: |a — &| < c < a + b. 
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(b) The evaluation (r,7)^ of a spin network (F, 7) is defined to be zero if it is not admissible. An admissible 
spin network is evaluated by the following algorithm. 

• Use the cyclic ordering to thicken the vertices into disks and the edges into untwisted bands. 

• Replace the vertices and edges by the linear combinations of arcs as shown in Figure O 

• Finally the resulting linear combination of closed loops is evaluated by assigning the value (—2)" to 
a term containing n loops. 




Figure 2. The rules for evaluating a spin network. Replace vertices and edges by linear combinations of 
arcs and count the number of closed loops. 

For example, in Figurc[T]the O network and the Cube are admissible, the 5-sided drum and the 9j'-symbol 
are not. Note that the admissibility condition is equivalent to saying that the strands can be connected at 
each vertex as in Figure [2l Note also that cubic ribbon graphs F are allowed to have multiple edges, loops 
and may be disconnected with components that contain no vertices. In addition, F may be and non-planar 
(contrary to the requirement of many authors |We| IMoul IKL| ). as long as one fixes a cyclic ordering of the 
edges at each vertex. The latter condition is implicit in |Pel| . It turns out that a different cyclic ordering 
on a vertex of a spin network effects its evaluation by a sign; see Lemma |6 . 1 1 below . 

1.2. A fundamental problem. It is easy to see that if (F,7) is an admissible spin network and n is a 
natural number, then (F, 727) is admissible as well. A fundamental problem is to study the asymptotic 
behavior of the sequence of evaluations (F,n7)^ when n is large. This problem actually consists of three 
parts. 

Problem 1. Fix an admissible spin network (F,7). 

(a) Prove the existence of an asymptotic expansion of the sequence (F,7i7)^ when n is large. 

(b) Compute the asymptotic expansion of (F,7i7)^. 

(c) Identify the terms in the asymptotic expansion with geometric invariants of the spin network. 

Although any spin network can be evaluated in terms of 3j and 6j-symbols, a positive solution to parts 
(a), (b) or (c) of Problem [1] for 3j and 6j-symbols does not imply a solution to for arbitrary spin networks. 

Problem [1] is motivated by the belief that the quantum mechanics of particles with large spin will approx- 
imate the classical theory. There is a large physics literature for Part (b); see for example [VMK . On the 
other hand, the literature for Part (a) is relatively new and short, and to the best of our knowledge concerns 
only 3j and 6j-symbols with certain labellings. Roberts used geometric quantization techniques to prove the 
leading asymptotic behavior of 6j'-symbols in the so-called Euclidean case; see [Rbll lRb2] . 

In order to get an exponentially (and not factorially) growing asymptotic expansion we will rescale the 
evaluations of spin networks by factorials of linear forms on the labels of their edges as follows. 

Definition 1.2. We define the standard normalization of a spin network evaluation to be 



(1) (F,7) = ^(r,7>^ 

where X! is defined to be the product 
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where the colors of the edges emanating from vertex ti, and V^(r) is the set of vertices of F. 

There are many useful properties of this evaluation: 

• (0,7) is still an integer (see Theorem [5]), 

• the sequence (F, 717) is exponentially bounded (see Theorem [1]), 

• the spin generating series is a rational function (see Theorem [S] below). 

2. Existence of asymptotics 

2.1. Asymptotics via G-functions. Our first result is a complete solution of Part (a) of Problem [1] for 
all spin networks, using the theory of G-functions. To state it, we need to recall the notion of a sequence of 
Nilsson type [Nij . 

Definition 2.1. We say that a sequence (a„) is of Nilsson type if it has an asymptotic expansion of the form 

(3) E ^"""(l0g«)''5A,a,/3/iA,a,/5(l/n) 

where 

• the summation is over a finite set, 

• the growth rates A are algebraic numbers of equal magnitude, 

• the exponents a are rational and the nilpotency exponents (3 are natural numbers, 

• the Stokes constants S\^a.f3 are complex numbers, 

• the formal power series h\^a,f}{x) € -f4r[[a;]] are Gevrey-1 (i.e., the coefficient of x" is bounded by 
G"n! for some G > 0), 

• if is a number field generated by the coefficients of hx^a.pix) for all A, a, /3. 

The above definition uses the notion of an asymptotic expansion in the sense of Poincare, given in jOj 
Sec.l]. In the special case where there is only one growth rate, the notation 



a„~A"r."(logn)''^^ 



means that for every e N we have 

/ Af-l 
lim a„A-"n-"(logn)-^ — 



In Section [7] we will prove the following. 

Theorem 1. For every spin network (F,7), the sequence (F, 717) is of Nilsson type. 

Theorem[T]follows from the theory of G-functions by combining Theorems[2]and[3]below. G-functions were 
introduced by Siegel ([Si]) in relation to transcendence problems in number theory. Many of their arithmetic 
and algebraic properties were established by Andre in |Anj . G-functions appear naturally in Geometry (as 
Variations of Mixed Hodge Structures), in Arithmetic and most recently in Enumerative Combinatorics. For 
a detailed discussion, see |Ga4j and references therein. Let us recall the notion of a G-function. 

Definition 2.2. We say that series G{z) — X^^o ^ G-function if 

(a) the coefficients a„ are algebraic numbers, 

(b) there exists a constant G > so that for every n S N the absolute value of every conjugate of a„ is 
less than or equal to G", 

(c) the common denominator of oq, . . . , a„ is less than or equal to G", 

(d) G{z) is holonomic, i.e., it satisfies a linear differential equation with coefficients polynomials in z. 

The following theorem is proven in Section [T] The proof is elementary given the results of |Ga4j and 
the fact that an evaluation (F, 717) can be written as a multi-dimensional combinatorial sum of a balanced 
hypergeometric term. 
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Theorem 2. The generating function 

oo 

(4) Fr,^(z) = ^(r,n7)z" 

ri=0 

of an arbitrary spin network (r,7) is a G-function. 

The next result follows by a combination of results of Andre, Katz and Chudnovsky-Chudnovsky; see [An] 
and also [Ga41 Sec. 2]. 

Theorem 3. If G{z) = X)riLo G-function, then (a„) is a sequence of Nilsson type. 
Every G-function has two basic invariants: an analytic and an arithmetic one: 

• the radius of convergence at z = 0, which equals to |A|~^ where A is a growth rate of (a„) as in 
Definition!^ 

• a number field K corresponding to (a„) as in Definition 12.11 

Theorem [2] motivates the following invariants of a spin network: its spectral radius and its number field. 
Definition 2.3. Let (r,7) be a spin network. 

(a) The spectral radius pr.-y is defined to be the inverse of the radius of convergence of the generating 
function Fy^^{z) at z = 0. 

(b) The number field K^^^ is defined to be the number field associated to the G-function (jH). 

(c) The spectral radius pY and the number field Ky of a cubic ribbon graph is defined to be pr,2 and 
Kv,2 respectively, where 2 denotes the admissible coloring of all edges of F by 2. 

2.2. Computation of asymptotics. The proof of Theorem [1] is not constructive. We know of several 
different methods for computing the asymptotic expansions of spin network evaluations. We will list them 
here briefly and illustrate them in the following sections. 

(a) The method of Borel transform for one-dimensional sums of balanced hypergeometric terms. This 
method has been developed for one-dimensional sums in |CG1[ ICG2i IGa3| ) and gives exact for- 
mulas for the complete asymptotic expansion of such sums. The method ought to extend to multi- 
dimensional sums, but this still needs to be developed. Since 6j-symbols are given by one-dimensional 
sums of a balanced hypergeometric term (in all cases, Euclidean, Plane or Minkowskian) , this allows 
us to compute exactly the leading asymptotic behavior of 6j-symbols. 

(b) The method of G- functions. We mentioned already that the generating series FY^.y{z) of ([5]) is a 
G-function in a non-constructive way, but more is actually true. In TheoremlHlwe show that FY.-y{z) 
is the diagonal of an explicit rational function. Thus, the asymptotics of the Taylor coefficients 
(r,n7) of i^r,7(^) can be computed in principle via Hodge theory, along the lines of BK ; see also 
|PW1[ IPW2] . Unfortunately, Pemantle- Wilson's results require genericity assumptions that do not 
apply even in the case of the ^3,3 or the Cube spin network. 

(c) The method of Wilf-Zeilberger. The fundamental theorem of Wilf-Zeilberger implies that for an 
arbitrary spin network (r,7), the sequence (F, 717) is a multi-dimensional sum of a hypergeometric 
term, and thus, holonomic, i.e., it satisfies a recursion relation with coefficients polynomials in n, 
|WZ1] . The WZ method is constructive, and has been implemented in several platforms, |PWZ1 
IPaRii |WZ1| . WZ has been implemented for multidimensional sums, but in practice it works best for 
one or two dimensional sums. Using the WZ method, we compute the asymptotics of three simple 
examples, a Euclidean 6j'-symbol, a Plane one and a Minkowskian 6j'-symbol, in Section [TUl 

(d) The ansatz for the asymptotics of hypergeometric multi-dimensional sums. In |Ga2] the first author 
developed a constructive ansatz for locating the singularities of the generating series of balanced 
hypergeometric multi-dimensional sums. The conjecture has been confirmed for one-dimensional 
sums using the theory of Borel transform. For higher dimensional sums the conjecture requires a 
genericity assumption (finiteness of the solution set of the so-called Variational Equations) which 
does not seem to hold for spin networks with more than six edges. A proof of the ansatz in the 
higher dimensional case is expected to follow by developing further the theory of Borel transform. 

The practical limitation of these methods is discussed in Section [T2l 
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3. ASYMPTOTICS OF 6j-SYMBOLS 

3.1. The geometry of the tetrahedron and asymptotics of the 6j-symbol. In the special case of the 
tetrahedron graph we present a complete solution to Parts (b) and (c) of Problem [1] using the method of 
Borel transform. We can explicitly compute the asymptotic expansions for all possible colorings, even in the 
degenerate and non-physical cases. Moreover the quantities in the asymptotic expansion can be expressed 
as geometric properties of the dual tetrahedron graph interpreted as a metric polyhedron. 




Figure 3. The three realizations of the dual metric tetrahedron depending on the sign of the Cayley-Menger 
determinant det(C). Also depicted are the two exponential growth factors A± of the 6j-symbol. 

In this section (^>t^, 7) denotes the tetrahedral spin network colored by an admissible coloring 7. Consider 
the planar dual graph (also a tetrahedron) with the same labelling on the edges. Regarding the labels 
of as edge lengths one may ask whether (X[;'',7) can be realized as a metric Euclidean tetrahedron with 
these edge lengths. The admissibility of ( /t\,7) implies that all faces of the dual tetrahedron satisfy the 
triangle inequality. A well-known theorem of metric geometry ^1] implies that ( 7) can be realized in 
exactly one of three flat geometries 

(a) Euclidean 3-dimensional space 

(b) Minkowskian space R^'^ 

(c) Plane Euchdean M^. 

Moreover, the above trichotomy is decided by the sign of the Cayley-Menger determinant det(C) which is a 
degree six polynomial in the six edge labels; see Definition l9.ll Without loss of generality, we will assume that 
i /K i'y) non-degenerate in the sense that all faces of ' are two dimensional (i.e., the triangle inequalities 
are strict). In the degenerate case, the evaluation of 6-symbol is a ratio of factorials, whose asymptotics 
are easily obtained from Stirling's formula. Notice that ( ,7) can be degenerate only in the Plane and 
Minkowskian cases. 

3.2. Numerical invariants of tetrahedra. Every sequence of Nilsson type comes with its own numeri- 
cal invariants: its growth rates, exponents, nilpotency exponents. Stokes constants and number field; see 
Definition 12.11 Theorem [T] implies that 6j-symbols have such invariants. In this section we will define a 
set of numerical invariants of tetrahedra and in the next section we will identify those with the numerical 
invariants of the corresponding 6j-symbols. In this section, we will fix a nondegenerate tetrahedral spin 
network ( A:, , 7) . 
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Definition 3.1. (a) Let 9a be the exterior dihedral angle corresponding to edge a of the metric reahzation 
of 

(b) Consider the growth rates A± defined by: 

(5) A± =e±*^"''"t. 

(c) Consider the number field 



(6) K = Q(V~det(C)) 

The next Lemma, proven in Section 19.21 explains the top part of Figure [31 

Lemma 3.2. The growth rates A± have the following properties. 

(a) A+A_ = f 

(b) In the Euclidean case we have additionally that A_ = A_|_ and A+ ^ A_ and consequently |A_|_| — 
|A_| = 1, and A± ^ ±1. 

(c) In the Plane case, we have additionally that A-(_ = A_ and consequently A± e {—1, !}• 

(d) In the Minkowsian case, we have additionally that A± G M \ {0}. Exactly one of A± is inside and 
one is outside the unit circle. Let A denote the one inside the unit circle. 

Definition 3.3. We define the exponents a± by 

if 7 is Eulcidean or Minkowskian 



(7) (a+,Q;_) = 



^ 5^ if 7 is Plane. 



2 ' 2 
4 _5 
I, V 3' 3 

Definition 3.4. (a) If 7 is Euclidean, we define 



(b) If 7 is Minkowskian, we define 
(9) 5± = T 



(c) If 7 is Plane, we define 
(10) 5i 



V6ttVo\ 
r(4/3)(12A)i/3 



The definition of A and Bj (homogeneous polynomials in 7 of degree 2 and 3 respectively) in the Plane 
case can be found in Section [ 



3.3. Asymptotics of 6j-symbols. In this section we will express the asymptotics of the 6j symbols (using 
the minor change (•)^ in normalization explained in Definition 19. 4p in terms of the numerical invariants of 
the previous section. 

Theorem 4. Fix a nondegenerate tetrahedron { A\.,j). 

(a) 7/7 is Euclidean or Plane, there exist h± S Js^[[a;]] so that 



(11) (^,n7)^'-A>"+5+/i+ l^^^+Aln'^-S^h^ (^^^ 



(b) 7/7 is Minkowskian, let S be the Stokes constant corresponding to A. There exists h G -?ir[[x]] so that 
Moreover, K is the associated number field. 



8 



STAVROS GAROUFALIDIS AND ROLAND VAN DER VEEN 




Equation (fTT|) implies in the Euclidean case that 

which is the well-known Ponzano-Regge formula for the asymptotics of Euclidean 6j-synibols. In the Eu- 
clidean and Minkowskian cases these asymptotic formulas were conjectured by Ponzano-Regge [PRj . An 
elementary but non-rigorous proof of the Euclidean case was given in |Gu| . The Euclidean case was rigor- 
ously proven by |Rbl| but there was no proof for the Minkowksian case. The formula for the Plane case is 
new. 

Remark 3.5. Notice that exponents a± in the Plane case of Theorem |4] is not a half-integer but rather the 
third of an integer. This phenomenon (accompanied by the collision of the two points in the Plane case of 
Figure [3]) is natural from the point of view of local monodromy of G-functions, and is often overlooked in 
the physics literature. 

An illustration of Theorem 3] in all three geometries is given in Section [TUl using the WZ method. 

Remark 3.6. Using Theorem |4l we can compute explicitly the leading asymptotics of an infinite family of 
spin networks whose underlying graphs can be created from the tetrahedron graph by repeatedly replacing 
a vertex by a triangle. As follows from the recoupling theory discussed in Section [S] the corresponding spin 
network evaluation will change by multiplication by a certain 6j-symbol and a theta. We thus find the 
asymptotics for all these networks. In particular, this leads to spin networks with an arbitrary number of 
growth rates, and with associated number fields which are the composite of an arbitrary number of quadratic 
number fields. 

4. A GENERATING FUNCTION FOR SPIN NETWORK EVALUATIONS 

In this section we discuss the generating series Sr of a spin network evaluation for a fixed underlying 
graph r. As it turns out, S'r is a rational function. Among other things, this implies that the G-functions 
of Theorem [2] come from geometry and more precisely from a Variation of Mixed Hodge structures, much in 
the spirit of Bloch-Kreimer [BK| . 

The next definition introduces the generating series Sr of a fixed graph F with cyclic ordering. Sr is a 
formal power series in commuting variables Ze where e G E{T), the set of edges of F. 

Definition 4.1. The spin generating function Sr of a cubic ribbon graph F is defined to be 

(12) 5r = ^(F,7) n eQ[[ze,eei?(F)]] 

7 e<£E(r) 

where the sum ranges over all colorings 7 of the graph F. 

We will express the spin generating function of an arbitrary graph in terms of the curve polynomial that 
we will define now. A curve in a graph F is a 2-regular subgraph. By definition the empty set is also a curve. 
The set of all curves in F is called Cr- 

Definition 4.2. Given a cubic ribbon graph F and A" C Gr we define 

(a) a polynomial 

(13) Pr,x= ^ ex(c)n^e. 

where ex(c) = —1 (resp. 1) when cG X (resp. c ^ X). 

(b) We will call Ppj the curve polynomial of F. 

With the above notation, we define l{X) to be the minimal number of intersections when X is drawn on 
the thickened ribbon graph F and 

(14) - ^ E 

YcCr 



ASYMPTOTICS OF CLASSICAL SPIN NETWORKS 



9 



The next theorem proves the rationaUty of Sr- 



Theorem 5. For every cubic ribbon graph T we have 



(15) Yl axPr,x^Hze,e 



eE{r)]]nQ{z,,eeE{T)). 



XdCr 



Corollary 4.3. For every spin network (r,7), the evaluation (r,7) is an integer number. 

Corollary 4.4. (We] When T is planar with the counterclockwise orientation, then ax = unless X is 
empty, in which case a0 = 1 . It follows that 



recovering an earlier theorem by Westbury [We) . 

The next theorem shows that the generating series FY^-f{z) of ([3]) is a diagonal of a rational function. Let 
us recah the notion of a diagonal. 

Definition 4.5. Given a power series f{xi, . . . , Xr) G Q[[a;i, . . . , Xr]] and an exponent / = (ii, . . . , v) G N+, 
we define the /-diagonal of / by 



Consequently, the G- function Fr^.y{z) is the diagonal of a rational function, and thus it comes from 
geometry; see Lemma 17.41 below. 

We end this section with the following curious proposition. 

Proposition 4.6. If T is a 2-connected cubic ribbon graph, then the evaluations (r,7) for all 7 uniquely 
determine F as an abstract trivalent graph. Moreover, we only need to use the colors and 1. 



In Section [6] we discuss the recoupling method for evaluating spin networks which reduces the evaluation 
of an arbitrary spin network to the 3j and 6j-symbols. Using this method, and an explicit formula for the 
3j and 6j-symbols, we evaluate the Drum and the K^^^ graphs. 

In Section [7| we prove that the evaluation of a spin network is given by a multi-dimensional sum of a 
balanced hypergeometric term. This, together with the results of (Ga4j . implies Theorem [2] and Theorem [T] 

In Section [8] we apply techniques of Borcl transform (developed in [CGll ICG2( iGa 3 ) that give exact 
formulas for the leading asymptotic behavior of one-dimensional sums of hypergeometric terms. Since 6j- 
symbols are given by such sums, this allows us to compute exactly the leading asymptotic behavior of 
6j-symbols. 

Section [9] identifies the combinatorial leading asymptotic behavior with geometric quantities associated 
to the metric tetrahedron in all cases, Euclidean, Plane or Minkowskian. 

In Section [10] we present an alternative constructive method to compute the asymptotic behavior of 
evaluations of spin networks, namely the WZ method. 

In Section [Tl] we discuss in detail the chromatic evaluation of spin networks and obtain the integrality 
and rationality results of the spin generating series. 

5.1. Acknowledgment. The results were conceived in a workshop in Aarhus, Denmark, and presented in 
HaNoi, Vietnam and Strasbourg, France in the summer of 2007. The authors wish to thank the organizers 
for their hospitality. S.G. wishes to thank D. Zeilberger for many enlightening conversations. 
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where [x'^^]{f) denotes the coefficient of x™^ . . . x™'' in /. 

Theorem 6. For every spin network (F,7) we have 

(17) Fr,^ = A^Sf 



5. Plan of the paper 
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6. Evaluation of spin networks 
We start by recording some elementary facts about spin network evaluations. 

Lemma 6.1. Let (r,7) be a spin network. 

(a) Inserting a vertex colored (0, a, a) in the interior of an edge of T colored a does not change the 
evaluation of the spin network. 

(b) Changing the cyclic ordering at a vertex whose edges are colored a, b, c changes the evaluation by a 
sign (-l)("(«-i)+''(''-i)+c(c-i))/2 

Proof, (a) The chosen normalization introduces an extra factor 1/a! for the new vertex labeled (0, o, a), 
while it follows from the definition that one also inserts an extraneous summation over permutations in the 
pre-existing edge labeled a. Since X^o-gS SreS sign(cr)(Tsign(r)T = a! J2aes sign((T)cr the evaluation is 
unchanged. 

(b) Changing the cyclic order at a vertex with edge labels a,b,c has the following effect. The alternating 
sum at each of the adjacent edges is multiplied by the permutation that turns the arcs in the edge by 180 
degrees. This element has sign a{a — l)/2 in Sa- □ 

As a consequence of part (a) of the above lemma, an edge labeled in a spin network can be removed 
without affecting the evaluation. 

There is an alternative bracket normalization {T,j)^ of the evaluation of a spin network (r,7) which 
agrees with a specialization of the Jones polynomial and the Kauffman bracket. 

Definition 6.2. The bracket normalization of a spin network (r,7) is defined by 

(18) (r,7)^ = ^(r,7>^ 

where 

(19) n ^(^)! 

eeE(r) 

This normalization has the property that it coincides with the Kauffman bracket (Jones polynomial) of a 
quantum spin network evaluated at A = — 1 [KL| . However, (r,7)-^ is not necessarily an integer number, 
and the analogous generating series does not satisfy the rationality property of Theorem [51 

Proof, (of Proposition ^^ For convenience we work with the unnormalized evaluation (F, 7)^. A 2-connected 
graph is determined up to isomorphism by its 2-regular subgraphs, simply take a polygon for every 2-regular 
subgraph and identify equal edges [Di]. To decide whether a set X of edges forms a 2-regular subgraph, let 
7 be the coloring that assigns 1 to members of X and to the other edges. By definition of the spin network 
evaluation X is an n component 2-regular subgraph if and only if (F,7)^ = (~2)". □ 

6.1. Evaluation of spin networks by recoupling. In this subsection we describe a way of evaluating 
spin networks by recoupling. We will reduce the evaluation of spin networks to multi-dimensional sums of 
6j and 3j-symbols. The value of the 6j and 3j-symbols is given by the following lemma of |KL] and [We| . 
using our normalization. 

Lemma 6.3. (a) Let ( 1,7) denote a tetrahedron labeled and oriented as in Figure [1] or Figure [5] (left) 
with 7 — (a, 6, c, d, e, /). Its normalized evaluation is given by 

minSj / h \ 

^ \Si~ k,S2 - k,S3- k,k-Ti,k~T2,k-T3,k-T4j 

fc— max i i 

where, as usual 

a 



ai,a2, . . . ,arj ai\ . 
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denotes the multinomial coefficient when ai + . . . aj. ~ a, and Si are the half sums of the labels in the three 
quadrangular curves in the tetrahedron and Tj are the half sums of the thee edges emanating from a given 
vertex. In other words, the Si and Tj are given by 

(21) Si = ^{a + d + b + c) S2^^{a + d + e + f) S3 = ^{b + c + e + f) 

(22) Ti = i(a + 6 + e) T2 = ^{a + c + f) = ^{c + d + e) n = ^{b + d + f). 

(b) Let (8, 7) denote the 8 spin network of Figure [1] admissibly colored by 7 = (a, b, c). Then we have 

(23) 0,7) = (-1)^ [^^^^ + 1) {^^^ 

Note that the value of a n-colored unknot is by Lemmas 16. II and 16.31 equal to (— l)"(n + 1). 

Recoupling is a way to modify a spin network locally, while preserving its evaluation. This is done as 
in the topmost picture in Figure [H This formula is called the recoupling formula and follows from the 
recoupling formula in [KL,, using our conventions. The other two pictures in the figure show two important 
special cases of the recoupling formula, the bubble formula and the triangle formula. 




Figure 4. The recoupling formula (top), the bubble formula (left) and the triangle formula (right). The 
sum is over all k for which the network is admissible, and Sk,i is the Kronecker delta function. 

The bubble formula shown on the left of Figure [3] serves to remove all bigon faces. Likewise the triangle 
formula can be used to remove triangles. To see why the recoupling formula suffices to write any spin network 
as a multi sum of products of 6j'-symbols divided by thetas we argue as follows. Applying the recoupling 
formula to a curve in the graph reduces its length by one. Keep going until you get a multiple edge which 
can then be removed by the bubble formula. 

Although the triangle formula follows quickly from the bubble formula and the recoupling formula it is 
important enough to state on its own. This formula shows that the evaluation of the class of triangular 
networks is especially simple. The triangular networks, mentioned in Remark 13. 61 are the planar graphs that 
can be obtained from the tetrahedron by repeatedly replacing a vertex by a triangle. By the triangle formula 
any triangular network is simply a product of 6j-symbols divided by thetas. 

To illustrate how recoupling theory works, let us evaluate the regular s-sided drum and 3. Consider 
the s-sided drum network as shown in Equation (I24p (for s = 5) where every edge is colored by the integer 
n. In the figure we have left out most of the labels n for clarity. By convention unlabeled edges are colored 
by n. Performing the recoupling move on every inward pointing edge we transform the drum into a string 
of bubbles that is readily evaluated. 
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Observing that if n is odd the network is not admissible (and thus evaluates to zero), and denoting the 
tetrahedron and the theta with one edge colored by k and the others by n by iS'(7i, k) and 9{n, k) we conclude 
the following formula for the n-colored s-sided drum. 

Proposition 6.4. If n — 2N is even we have 

and if n is odd we have (Drums, n) = 0. 

For small values of s the drum can be evaluated in a more straight forward way, thus providing some 
well known identities of 6j'-symbols. Namely when s = 1 we get zero, when s = 2 we find some thetas and 
when s = 3 we have by the triangle formula a product of two 6j-symbols thus giving a special case of the 
Biedenharn-EUiott identity KL]. For s — A we find a formula for the regular cube, we know of no easier 
evaluation in this case. 

A similar computation for K^^s cyclically ordered as a plane hexagon with its three diagonals gives the 
following. 

Proposition 6.5. If n ^ 2N is even we have 




and if n is odd we have {Kj,,^,n) = 0. 

Note the similarity between Drums and The only difference is the sign that comes up in the 

calculation when one needs to change the cyclic order. The extra sign makes (1^3,3, 2N) — for all odd N. 
This is because changing the cyclic ordering at a vertex takes the graph into itself, while it produces a sign 
(— 1)^ when all edges are colored 2N. 

7. G-FUNCTIONS AND EXISTENCE OF ASYMPTOTIC EXPANSIONS 

7.1. Existence of asymptotic expansions. In this section we prove the existence of an asymptotic expan- 
sion for all spin networks using the theory of G-functions and the combinatorial fact that any spin network 
can be evaluated as a multi-dimensional sum of a balanced hypergeometric term. 

Definition 7.1. An r-dimensional balanced hypergeometric term tn,k (m short, balanced term) in variables 
(n, fc), where n e N and k — (fci, kr) (z , is an expression of the form 

r J 

(25) t,,^k = C-l[CJ^l[Mn,ky^^ 

i=i j=i 

where Ci are algebraic numbers for « = 0, . . . , r, — ±1 for j — 1, . . . , J, and Aj are integral affine linear 
forms in (n, fc) that satisfy the balance condition 

J 

(26) ^ ejAj — constant 
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We assume that the convex polytope Pi in W defined hy Aj{l,w) > for j — 1, . . . , J and w = {wi , . . . , Wr) 
is compact. 

Lemma 7.2. For every spin network (r,7) there exists a balanced term t{n,k) with k ~ (fci, . . . ,kr) such 
that 

(27) (r,7i7)= ^("'^) 

for every n G N. 

Proof. Using the recouphng formulae from the previous paragraph we can write (F, 7) as a multi-dimensional 
sum of products of 6_;/-symbols, 3_;/-symbols and unknots (i.e., Ij-symbols) with a denominator consisting of 
3j-symbols. It follows from Equations ^U\j and that the 6j-symbols (resp. 3j-symbols) are balanced 
1-dimensional (resp. 0-dimensional) sums, thus the ratio of the product of the 3j-symbols by the product of 
the 3j-symbols is a balanced multi-dimensional sum. The unknots can be written as {—l)'^{k + l)l/kl and 
are therefore balanced as well. It is easy to check that admissibility guarantees that the multi-dimensional 
sum has finite range. □ 



Beware that the term t(n, k) constructed in the proof of lemma [7T2] is neither unique nor canonical in any 
sense. Lemma [7.21 and the following result implies Theorem [TJ 

Theorem 7. |Ga4t Thm.2] (a) If (a„) is a multi-dimensional sum of a balanced term, then the generating 
series 



F{z) = ^ QnZ 



is a G-function. 

(b) If F{z) is a G-function, its Taylor series at z = is a sequence of Nilsson type. 

Remark 7.3. This argument also implies that the sequences (r,?i7)'^ and (r,n7)^ are also asymptotic 
to a rational power of n\ times a sequence of Nilsson type. This follows from the fact that the above 
normalizations of the spin networks differ from Definition 11.21 by a ratio of factorials of linear forms in n, 
which can be asymptotically expanded using Stirling's formula. 

Let us end this section with a lemma which shows that diagonals of rational functions are G-functions that 
come from geometry. This is relevant to spin networks since we know from Theorem [5] that the generating 
function i^r,7 for the evaluations of a spin network is a diagonal of a rational function. 

Fix a power series f{xi, . . . ,Xr) G Q[[a;i, . . . , a;^]] convergent at the origin and an exponent I ~ (ii, . . . , v) G 
N!j_ and consider the diagonal {Ajf){z) G Q[[z]] as in Definition 14.51 Let C denote a small r-dimensional 
torus around the origin. Then we have the following. 

Lemma 7.4. With the above assumptions, we have 

(28) (A,/)(z) ^ -i- / li^lll^^dx, A • • • A dxr. 

[ZTny Jc xi . . . Xr — z 

Proof. With the notation of Definition llGl an application of Cauchy's theorem gives for every natural number 
n 

_ 1 [ f{xi,...,Xr) 



Summing up for n and interchanging summation and integration concludes the proof. □ 

If in addition /(xi, . . . , Xr) is a rational function, then the singularities of the analytic continuation of the 
right hand-side of ((28l) can be analyzed by deforming the integration cycle C and studying the corresponding 
variation of Mixed Hodge Structure as in [BKj . Such G-functions come from geometry; see [Ani IBK) . 
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8. Computation of asymptotic expansions using Borel transform 

The results of Section [7] are powerful but non-constructive. In this section we will use constructive results 
from the theory of Borel transform, known in the case of balanced terms of rank 1. These results appear 
in [CG1|, ICG21 IGa4| . and in the present paper we will quote them without proofs. Let us specialize to 
1-dimensional balanced terms tn,k, that is expressions of the form 



J 

(29) tn,k^CSC'^l[A,{n,k)l'^ 

satisfying the balancing condition (l26l) . The relevant example of these is the 6j-symbol, Equation (|20p . 
Fixing a 1-dimensional balanced term summation over k (it is a finite sum) gives rise to a sequence 

(30) a„ = y^JnM 

fcez 

To describe and compute the asymptotics of such sequences we need the following notation. For afhne linear 
forms we use the notation A(n, k) — VQ{A)n + vi{A)k -t- v{A) and A{u) = wo(^) + vi{A)u. Next we define 
the multivalued analytic functions 

J 

(31) V{u) = logCo+7/logCi+^e,A,(7.)log(A,(u)) 

J 

(32) W{u) = [|Aj(u)'^/2+.,t,(A,) 

V and W are analytic functions on the open interval (m, A'l) (which is determined by Aj, j = 1, . . . , J) and 
they have single-valued analytic continuation in the doubly cut complex plane C\ ((— oo,m] U [M,oo)). V 
and W can be further analytically continued as multivalued functions in the complex plane minus a finite 
set of points C \ {""nj Aj{u) = 0}. 
Define the constants 



(33) f^=\^^i' ^ = ^<^AAj), 

i=i i=i 

The set of critical points of V satisfies the Variational Equation 

J 

(34) Ci]^ylj(u)^^-^i(^^) = 1 

i=i 

If u is a critical point of V , then e^(") — \{u) where 



(35) A(u)= n ^^(^ 



j:Aj(«)#0 



Let K denote the number field generated by the solutions of the Variational Equation. 

Suppose that u is a critical point of order m for V ^ that is V'^^^ (u) — for j — 1, . . . ,m and (u) =/= 

and W is non-vanishing at u, i.e., W{u) ^ 0. Let 

(36) c™ = 2rf^^Vm-|-l)!^. 
and let 

(37) a{u) = ^ + 1/ + 1 5— - 

m + 1 



ASYMPTOTICS OF CLASSICAL SPIN NETWORKS 



15 



and 

(38) Sin) = (27r)^c„— — ^^^1-^ 

We will require a modification of the above definition of the Stokes constant when u lies in one of the cuts 
(— oo, m] U [M, oo). In that case, we define 

(39) S{u) = (2^)^c„ l^^""^' 



(|l/(™+i)(u)|)" 
Theorem 8. With the above assumptions we have 

(40) an^Y. KuTn"-S{u)hu (- 



where G -^^[[2;]] is Gevrey-1 and u runs through a subset of solutions of the Variational Equation. The 
selected subset can be described algorithmically and in case of the 6j-symbols it is discussed below. 



9. ASYMPTOTICS OF THE 6j-SYMBOL 

The purpose of this section is to prove and explain Theorem 2] on the explicit expansions of the 6j- 
symbols. Since the asymptotics are phrased in terms of geometric properties of the tetrahedron dual to the 
spin network, we start with a section on the geometry of the tetrahedron. Next we apply Theorem [5] to find 
the asymptotic expansion and finally we work out the details in each of the geometric cases. 



9.1. The geometry of the tetrahedron. The asymptotics of the 6j-symbols turn out to be intimately 
related to the geometry of the planar dual tetrahedron. From now on, we will assume that 7 is non-degenerate, 
i.e. that the faces of (^^,7) are two dimensional triangles. 

Recall that the 6j-symbol is a tetrahedral spin network ( 7) admissibly labeled as in Figure [1] with 
7 = (a, 6, c, e, /). Its dual tetrahedron is also labeled by 7. The tetrahedron and its planar dual, 

together with an ordering of the vertices and a coloring of the edges of the dual is depicted in Figure O 
When a more systematic notation for the edge labels is needed we will denote them by dij (or ij if there 
cannot arise any confusion) it follows that 

(a, b, c, d, e, /) = (^12,^23,^14, ^34, ^13,^24) 




Figure 5. The tetrahedral spin network and its planar dual with a vertex and edge labeling. 

Recall from Section 13.1! that we interpret the labels of the dual tetrahedron as edge lengths in a 
suitable flat geometry. One may now ask for a condition that allows one to realize ( 7) in a flat metric 
space such that the edge lengths equal the edge labels. Labeling the vertices of ( 7) as in FigurelHlwe can 
formulate such a condition in terms of the Cayley-Menger determinant. This is a homogeneous polynomial 
of degree 3 in the six variables a^, . . . , f^. We give a definition of the determinant following [Kol| . 
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Definition 9.1. Given numbers dij define the (modified) Cayley-Menger matrix (C^) by C.y = 1 — 

for i, J > 1 and Cy = sgn(i — j') for when i = or j = 0. In terms of the coloring 7 = (a, 5, c, c?, e, /) of a 

tetrahedron, we have 



(41) C = 



We can now state the criterion for the realizability of (\[ ^7). For a proof see [Bl] or jKolj . 

Proposition 9.2. Consider a non- degenerate tetrahedron 7)- The sign of the Cayley-Menger deter- 
minant determines in what space the tetrahedron can be realized such that the edge labels equal the edge 
lengths 

• //det(C) > then the tetrahedron is realized in Euclidean space M.^. 

• J/det(C) = then the tetrahedron is realized in the Euclidean plane R^, 

• //det(C) < then the tetrahedron is realized in Minkowski space R2>i_ 

In each case the volume of the tetrahedron is given by Vol = ^ -y/j det(C)|. 

The 6 dimensional space of non-degenerate tetrahedra thus consists of regions of Minkowskian and regions 
of Euclidean tetrahedra. It turns out to be a cone that is made up from one connected component of three 
dimensional Euclidean tetrahedra and two connected components of Minkowskian tetrahedra. The three 
dimensional Euclidean and Minkowskian tetrahedra are separated by Plane tetrahedra. The Plane tetrahedra 
also form two connected components, representatives of which are depicted in Figure[3l The tetrahedra in the 
Plane component that look like a triangle with an interior point are called triangular and the Plane tetrahedra 
from that other component that look like a quadrangle together with its diagonals are called quadrangular. 
We use the same names for the corresponding Minkowskian components. An integer representative of the 
triangular Plane tetrahedra is not easy to find as the smallest example is (37, 37, 13, 13, 24, 30). 

We finish this section with a short discussion of the dihedral angles of a tetrahedron realized in either of 
the three above spaces. The cosine and sine of these angles can be expressed in terms of certain minors of the 
Cayley-Menger matrix. Define the adjugate matrix ad(C) whose ij entry is (—1)'+^ times the determinant 
of the matrix obtained from C by deleting the i-th row and the j-th column. Define Cujj to be the matrix 
obtained from C by deleting both the i-th row and column and the j-th row and column. 

The Law of Sines and the Law of Cosines are well-known formulae for a triangle in the Euclidean plane. 
The next Lemma is the Law of Sines and the Law of Cosines for a tetrahedron in all three flat geometries. 

Lemma 9.3. Let 9ki be the exterior dihedral angle at the opposite edge ij. The following formula is valid 
for all non- degenerate tetrahedra 



.. ad{C)ki + det(C) det{Ckkii) 

exp(Okii) — , 

Vad(C)fefead(C)H 

Proof. See |Kol| for a proof in the Euclidean and Plane case. The Minkowskian case is entirely analogous 
once one uses the right conventions on angles |SP| . Given a Minkowskian tetrahedron the angle between two 
outward pointing normal vectors has the form m: -\- iuj, where n G N and w G M. Since the normal vectors 
are time-like they can either be in the same sheet or in different sheets of the unit hyperboloid. If they are 
in the same sheet then oj is defined to be the length the hyperbolic line between them and n = 0. If they 
are in opposite sheets then —uj is the length of the line between one and the reflection of the other in the 
horizontal plane and n = 1. □ 

9.2. Asymptotics of the 6j-symboL We now start to apply Theorem [5] to determine the asymptotics of 
the non-degenerate 6j-symbol. In this section it is convenient to use the following normalization (r,7)^ for 
spin networks. 
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Definition 9.4. We define the unitary normalization (T,^) of a spin network evaiuation (r,7) to be 
wfiere 

9(7)= n V\{Q,<^v,b,,c.)\ 
vev{r) 

and ay, by, Cy are tfie colors of the edges at vertex v. 

As it turns out, with this normalization the asymptotics of the 6j-syinbols matches the corresponding 
geometry well. Since the theta is a quotient of factorials one can use Stirling's formula [O] to find its 
asymptotics and hence one can see how the new choice in normalization effects the asymptotic expansion. 
For later use we record the asymptotic expansion of the normalization factor in the case of a 6j-symbol. 

Let us define 

_ i-a + b + c)(-°+''+^)/-^(q -b + c)(°-^+^)/-^(q + b - c)(°+''-^)/4 
" (a + 6 + c)("+''+-)/4 

_ ((-g + b + c){a-b + c){a + b - c)(a + b + c))^'^ 
a + b + c 

Lemma 9.5. Let 7 be an admissible coloring of a tetrahedron. Then we have 



6(717) ^ ^ {2TT)'^SabeSacfScedSbdf{^abe^acf^ced^bdf)"' ( ^ ^ 



where fiQ — I and fii G Q(7) for all I. 

Thus, the asymptotics of (/;l\,7)'^, is reduced to the study the asymptotics of Since this is a 

1-dimensional sum of a balanced term given by Equation (P(7| , we can apply Theorem|S]to find its asymptotic 
expansion. Equation (|20p implies that the balanced term is given by 

^ ^ "''"n.(-^.-fc)!n,(fc--7^.)! 

With the notation of Theorem [H] we have J = 8, and the linear forms Aj{n, k) are given by 

Ai{n, fc) = fc + 1, A2{n, k) ~ nS\ — k, A^ln, k) = nS2 — k, A4^{n, k) ~ nS^ — k, 
^1.5(71, k) — k — nTi, Ag(n, k) — k — nT2, ^7(71, k) — k — nT^, As{n, k) = k — nT^ 
and ei = 1, €2 = . . . £8 = — 1- It follows that 

(43) = 1, /I = -3. 
The functions V and W from Theorem [S] in this case are 

(44) V{u) =ulog{-l) +ulogu- Y,iS^ - u) log(5, ~ u) - ^(li - T,) log(u - T,) 

(45) Win) ^ u'/'illiS, ~ u) l[iu - T,))''/^. 
V and W are analytic functions on the open interval 

(m,M) = (max{0,rj},min{5J) 

and they have single- valued analytic continuation in the doubly cut complex plane C \ ((—00, max{0, Tj}] U 
[minjS'i}, 00)). V and W can be further analytically continued as multivalued functions in C \ {0, Si,Tj}. 
The possible critical values of V{u) that determine the growth rate are given by 

(46) A, = 
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here u is a solution of the variational equation which is the equation 

(47) -uYl{S,-u)l[iu^T,)-' = 1 

At first sight, the above equation appears to be a quartic equation for u. However, the coefficients of u** 
and is zero because '^Sj — ^Ti. Thus, Equation (|T7l) is quadratic in u, with coefficients which are 
polynomial in a, . . . , /. It is convenient to introduce a fourth number 6*4 = 0. We can then write the 
Variational Equation P7)) in the form 

(48) E{u) = - JJiSi -u) + Y[{u - T.j) =Au^ + Bu + C' = Q 

where A, B and C" are homogeneous polynomials in 7 of degrees 2, 3 and 4 respectively. A calculation gives 
that 

(49) A = i(arf + 6c + e/) 

(50) B = -^{bc{b + c) + ad{a + d)+ef{e + f) + abc + abd + acd + bcd + abe 

+bce + ade + cde + acf + bcf + adf + bdf + aef + bef + cef + def) 

The next Lemma identifies the number field generated by the Variational Equation ()48|) in terms of the 
Cayley-Menger determinant. 



Lemma 9.6. Let D be the discriminant of the Variational Equation (j47p . Then we have 
(51) Z? = -det(C)/4 

where C is the Cayley-Menger matrix. The number field K generated by the solutions to the Variational 
Equation (|48p is given by 



(52) K = g(V-det(C)). 

Proof. The proof follows from an explicit calculation, keeping in mind that D and det(C) are homogeneous 
polynomials in 7 of degree 6. □ 



The next Lemma determines the order of vanishing of V at the solutions of Equation (|47p. 

Lemma 9.7. Let D be the discriminant of the Variational Equation (I47p . 

(a) The tetrahedron is degenerate iff there is a solution to the variational equation that equals Si, or Tj 
for some i, j . 

(Jo) Assume the tetrahedron is non-degenerate. The two solutions to the variational equation 

satisfy 

(c) // in addition D — then the two solutions coincide, V"{v±) — and we have 

l[iv~T,)V"'{v)^^2A 



where A is given in 

(d) In the Euclidean or Plane case, v± does not lie in the double cut (— cx3,to] U [AI,oo). In the 
Minkowskian case, v± are real and exactly one of them (call it v) lies in the double cut (— cx),m] U 
[M,oo). 

Proof. In the proof it is convenient to define 54 = so that all indices run from 1 to 4. This convention 
will only be used in this proof, (a) If the tetrahedron is degenerate then there are i,j such that Si = Tj. 
Plugging in this value for u will set both terms in the equation to zero as required. On the other hand if 
V = Tj is a solution to the variational equation then J^j(>5'i ^Tj) = 0, so that Tj = Si for some i and likewise 
when V = Si. 

(b) Note that A = (ad + bc + ef ) /2 > since all edge labels are positive so that a solution always exists. 
By (a) we know that the solutions to the variational equation are not equal to Si or Tj. 
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By differentiating V we find that —V"{v±) — J2 s -v± ~^ v±-t- ■ Using the variational equation 
Yii^i ^ u) = IK"*^ ~ ^j) '^^^ write 

(53) - n(-± - T.)v"i-±} = Uis^ - v±) E + n(-± - ^0 E 

We see that - n(^± ^ Tj)V"{v±) = E'{v±). Now choose the solutions such that v± = {-B ± ^/D)/2A so 
that E'{v±) = 2Av± + B ^ ±/D. 

(c) Since D — the solutions v± coincide and will be written as v. Differentiating we find V"'{v) = 
— Yli^i — + ~ the other hand D = and we assumed that our tetrahedron is non- 

degenerate, so by part (a) V"{v) — 0. Therefore 

1 I + y I = X>)=o 



{Si-v)^ ' j^. {St - v){Sj - v) ' {Si - v){v - Tk) St - V 
Adding these to similar equations for implies 

V'" («) = 2 y 

Multiplying through by the variational equation Y\{v — Tj) = Yli^i — v) we get 

n(^^ - T,)V"'iv) = 2E(^. - v)iSj -v)-{v- T,)iv - T,) = -2A, 

where A is given in Equation (j49|) . 

(d) The Euclidean or Plane case follows from (b). For the Minkowskian case, there are connected compo- 
nents (types) of Minkowskian tetrahedra: the rectangular type where Tj < v± and the triangular type where 
we have Ti,T2,r3 <v± < T4 where T4 is the distinguished triangle of the triangular tetrahedron. This fact 
is in turn proven by continuity. It is true by explicit computation (using the WZ method of Section [T0|) for 
the two examples 

(3, 4, 4, 3, 5, 5), (37, 37, 13, 13, 24, 30) 

of rectangular type and triangular type respectively. The result follows by continuity for all Minkowskian 
nondegenerate tetrahedra, since v± = Tj only for degenerate tetrahedra. □ 



From now on we will always assume that our 6j-symbols are non-degenerate. Using Lemma 19.51 and 
Theorem [8l it follows that the possible growth rates in the asymptotic expansion of ( ]; , , 7) ^ are given by 

A± = )^abe^acf^ced^bdf^v± 

We close this subsection with a lemma connecting these growth rates to the dihedral angles in the tetra- 
hedron. 

Lemma 9.8. Let 9a be the dihedral angle corresponding to the edge labeled a as introduced in Section \9.1\ 

(a) The growth rate A± satisfies 

(54) A±=e±*^»^"* 

(b) Moreover, 



(55) {SabeSacfScdeSbdfWt{v±)f 



Y{,{v±-T,) 

Proof, (a) The strategy is to collect equal powers in A±. It follows that we can write A± as a product 



i<3 

where we define 
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(56) 

^ _ {v± - (»j + tk + jk)/2){v± - {ij + il + jl)/2){ij - tk + jk){ij + tk-jk){ij ~ + + il - jl) 

~ {{ij + kl + ik + jl)/2 - v±){{ij + kl + jk + il)/2 - v±)VKlK^ 

and 

Ki — + ik + jk){ij — ik + jk){ij + ik — jk){ij + ik + jk) 
Now we will show that hij± = e^^''''\ Since det{Ckkii) — dfj and 4adH(C) — Ki we can clear denominators 
in the expression for e^^*"'* from Lemma 19.31 and the expression for hij±. Equality of the numerators is then 
shown by direct calculation using Equation (jSTjl . Equation ((54)) follows. 



(b) To show Equation (|55|). first note that sij ik jk — ij^ll^jj^ and hence 



{SabeSacfScdeSbdfyyt{v±)) = 



n(^. - «±) n(«± - mr.T^T.T.r 

In the product of the hij we can also collect like terms together in the form Kj in the denominator. We 

get 

ii ^3 {T^T2T:inYK^K^K:iKiY{{S,~v±Y 
Using the variational equation and hij± = g^s^ii proves part (b). □ 



Lemmas 19.81 and 19.71 prove Lemma 13.21 regarding the properties of the growth rates. 

Combining our results, we can finally prove Theorem |4] by reading off the numerical invariants from 
Theorem El First of all the growth rates A± are identified with the dihedral angles by Lemma 19.81 In the 
Euclidean and the Plane case both rates contribute, whereas in the Minkowskian case we include only the 
contribution of A, corresponding to the distinguished solution v of the Variational Equation that lies in the 
cut; see part (d) of Lemma [9.71 In the terminology of Section [5TT] and Lemma [9.71 we have A = A+ in the 
triangular component and A = A_ in the quadrangular component, and always |A| < 1. 

Next the exponent in the leading power of n is given by Equation ()37p and the order of vanishing is 
determined by Lemma l977l We have m = 1 in the Euclidean and Minkowskian cases and m = 2 in the Plane 
cases. 

Lemma |9 . 61 identifies the number field generated by the solutions to the Variational Equation with K from 
Definition 13.11 

Finally we consider the Stokes constants in all three cases. The Stokes constant is computed from Theorem 
[8] using Lemmas 19.81 and 19.71 In the Euclidean case we have to refine Lemma 19.81 (b) to choose the right 
square root of the left hand side. This is done by considering the simple example of the regular tetrahedron 
and extending by continuity. We thus find 



yw>i) v-mv±-T,)v"{v±) 

Now applying part (b) of Lemma 19.71 and Lemma 19.61 we can rewrite the above as 
S± = (27r)-i/2^!f _ £ 



In the Minkowskian case we have to choose the Stokes constant corresponding to A. This is computed as 
in the Euclidean case except that the solutions to the variational equation lie in the cut. We therefore use 
Equation ([39]) so that we take absolute values everywhere. This explains the sign in (l9|). It makes the Stokes 
constant positive by canceling the real part of the Minkowskian angle sum, which is Sir in the triangular case 
and 47r in the quadrangular case. 

Finally in the Plane case we have by Lemma [9771 that v ^ Tj and we know that V"{v) = 0. Furthermore 
C2 = r(4/3)6i/3, so 
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£(4/3)6^ \Wiv)\ _ r(4/3)6i/3|e'5:.e„/2| ^ r(4/3)6i/3 



\v"'{v)\-s ^in,(^^"r,)|i/2|y"'(„)|^ ^in,(^'-r,)|i/6(2^)i/3 

where we used part (c) of Lemma l97fl and the fact that the dihedral angles are real. Since the tetrahedron is 



Plane, we have y/D = 0, so u = —B/2A with the notation of (|48|) . Therefore we can take out the 2 A terms 
to get 

r(4/3)(12A)V3 
7rn|S + 2ylTj|i/6 

Recall from Equations ((22| . (|49| and ([50]) that Tj, A and i? are homogeneous polynomials of 7. A calculation 
shows that Bj := B + 2ATj is a homogeneous polynomial of degree 3. It is unclear to us what is the geometric 
interpretation of the rational function ^^/Jlj^i ^ Q(7)- 

10. Computation of asymptotic expansions using the WZ method 

In this section we will compute the asymptotic expansion of three representative 6j-symbols using the 
WZ method. The main observation is that a 6j-symbol is given by one-dimensional sum of a balanced 
hypergeometric term; see Equation ([20]) and Definition 17.11 The fundamental theorem of Wilf-Zeilberger 
implies that such sums are holonomic, i.e., they satisfy linear recursion relations with coefficients polynomials 
in n. The WZ method is constructive and works well for onc-dimensional sums, implemented in several 
platforms, see }PaRi| [PWZ| IWZ"T] . Given a recursion relation for the evaluation of the 6j-symbols, together 
with the fact that it is a sequence of Nilsson type, allows us to compute explicitly its asymptotic expansion. 

The three examples we consider are the simplest representatives of the three geometric types of 6j-symbols. 
Let 

7Euclidoan = (2, 2, 2, 2, 2, 2), 7Plano = (3, 4, 4, 3, 5, 5), 7Minkowskian = (4, 4, 4, 4, 6, 6). 

and consider the sequences (a„), (6„) and (cn) defined by 

n!6 ^ (~l)'=(fc + l)! 



A'^TEuclidean)'' - .3^ ^ ^M2 J2 (fc 



(3n + 1)!2 ^^^^ {k - 3n)!'*(4n - k) 



13 



^ nP(2n)!^(3n)!^ ^ (-l)'=(fc + l)! 



(6n + 1)!2 {k - 6n)!4(7n - fc)!(8n - fc)!(9n - A:)! 



(An. ^!^(3^)!^ V (-l)^(fc + l)! 

^ii^,n7Minkow«kia„; - ^ ^^,3 _ 7n)!4(8n - /fc)!(10n - /fc)!2 



The command 



<< zb.m 

loads the package of |PaRi j into Mathematica. The command 

teucl[n_, k_] :=n'.''6/ (3n + l) (-Ij-^k (k + 1) ! 1/ ((4n-k) (k-3n) !''4) 

defines the summand of the sequence (a„), and the command 

Zb[teucl[n, k] , {k, 3 n, 4n}, n, 2] 

computes the following second order linear recursion relation for the sequence (a„) 

-9(l + n) (2 + 3n)^ (4 + 3n)^ (451 + 460n + 115n2)a[n] + 

(3 + 2n) (319212 + 1427658n + 2578232 n2 + 2423109 + 1255139 n^ + 340515 n=+37835 n'') a[l+n] - 
9 (2 + n) {5 + 3n)^ (7 + 3n)^ (l06 + 230n + 115n^) a[2 + n] ==0 

This linear recursion has two formal power series solutions of the form 
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1 ^„ / -432±31iV2 109847 T 22320iV2 -18649008 ± 4914305i\^ 
°^±." - 576^^^ ^ 331776n2 ^ 573308928n3 

14721750481 ± 45578388960iv^ -83614134803760 ± 7532932167923iv^ 
^ 660451885056n4 38042028579225^^ 

-31784729861796581 T 212040612888146640^^2 
^ 657366253849018368n6 




where 

329 T 460i\/2 

A+ = 

729 

are two complex numbers of absolute value 1. These growth rates A± computed with the WZ method match 
those of Theorem ID Indeed, the Theorem states that 



g^i6 arccos(l/3) 

An algebraic calculation confirms that 

329±460iV5 



^±i6 arccos(l/3) 



729 



The coefficients of the formal power series aj-.„ are in the number field K — Q{\/ —2) and the Cayley-Menger 
determinant is det(C) = 2^ confirming Lemma WM and Theorem [H 

More is actually true. Namely, the sequence (a„) generates two new sequences „) and (/^_ „) defined 

by 



1 ' — ' 



1=0 



where /i±,o — 1- Each of the sequences (/i±,n) are factorially divergent. However, the generating series 
X^J^o -^"^i-n+i/"-' ^''^ G-functions (as follows from ^^id the sequences (/x±^„+i/n!) are of Nilsson 

type, with exponential growth rates K± — Aip. The asymptotics of each sequence gives rise 

to finitely many new sequences, and so on. All those sequences span a finite dimensional vector space, 
canonically attached to the sequence (a„). This is an instance of resurgence, and is explained in detail in 
pMl Sec.4]. 

The second order recursion relation for the Plane and the Minkowskian examples has lengthy coefficients, 
and leads to the following sequences (6±,n) and (c±^„) 



ASYMPTOTICS OF CLASSICAL SPIN NETWORKS 



23 



1 / 1 3713 25427 9063361 109895165 

b+,n = -TT^K 1 - — 



4/3 +\ 'Sn 46656^2 2239488^3 17414258688^4 lG4485552128n5 
1927530983327 



2437438960041984?i6 



1 / 37 3883 13129 5700973 14855978561 

-A" 1 - 



n5/3 - \^ 96n 46656n2 4478976n3 8707129344n4 3343537668096^5 
2862335448661 



+ 



2437438960041984n6 



1 ^„ / 336t 1369\/2 1769489 t 831792^2 67925105712^66827896993^2 
^±." - \ ^ 4032n ^ 1806336n2 ^ 21849440256^3 

5075437500833257 T 2589265090380768^2 

176193886224384^4 
100978405759997442992 t 98904713360431641651v^ 

552544027199668224n5 
685103512739058526758457 T 349782631602887151717776V5 \ 



+ 



247539724185451364352716 ' I 

where in the Plane case we have 

A_=A+ = -l, K = Q, det(C)=0 
and in the Minkowskian case we have 

696321931873- 111529584108\^ 

A+ = ^ = 0.794127... 

+ 678223072849 

696321931873+ 111529584108v^ 

A_ = ^ = 1.25924... 

678223072849 

K ^ QiV2), det(C) = -2^3^ 
again in agreement with Theorem [3) 

11. Chromatic evaluation of spin networks 

The purpose of this section is to prove the rationahty statement of the spin generating function stated in 
Theorem [S] of Section H) We will use the chromatic evaluation method which goes back to |Pe2j . Our proof 
builds on earlier work by |We] and |KL] on planar spin networks. 

The chromatic method for the evaluation of a spin network is based on introducing a generalized evaluation 
(r, 7)^ following the algorithm of Definition II. 1[ with the exception that a collection of n loops is assigned 
the value iV" instead of (—2)". Therefore, 

(r,7>^2 = (r,7>^. 

Now the crucial observation is that (r,7)^ is a polynomial in N with integer coefficients, so that (r,7)^ is 
determined once we know the values of (r,7)^ for all positive N. 

For positive N there exists an alternative interpretation of (r,7)^ as a contraction of tensors. For 
convenience we will orient the edges of T. Let y be an TV dimensional vector space and for a d Sn define 
the following permutation operators cr : F®" y^m^ a{ei^ (8) ■ ■ ■ (8) e^^ ) = e„(^i^-^ ■ ■ ■ <S) Sa-(i„)- At every edge 
colored n we can now speak of the map X^ctgs ^S^i'^)'^- Note that the single arcs in the expanded diagram 
for (r, 7) now correspond to copies of V. The systems of arcs at the vertices indicate how to contract 
(compose) the tensors on the edges. Since all the tensors involved are linear combinations of Kronecker 
deltas we can for each choice of permutations count the number of contracted deltas and this will exactly 
contribute a term ztiV*^ for some k since the trace of a delta is simply N. This coincides with the counting 
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and evaluation of loops in every state, since a contracted delta represents a loop. However there is also 
another way to carry out this tensor calculation. 

We can also calculate (r,7)^ by writing out the contraction of tensors into the single contributions of 
basis elements in the various tensor products. This will yield an explicit formula for the evaluation of a 
general spin network in terms of so called curve configurations on T. Recall that a curve is a 2-regular 
subgraph of F and a curve configuration is a function from the set C of all curves to the natural numbers, 
C : C ^ N. For such functions we define |£| to be the sum of the function values. Define sign (— l)'^'''^ 
where for any subset l{C) is the minimal number of intersections one can have when drawing C{c) copies of 
c for every curve on the thickened ribbon graph F. 

Lemma 11.1. For every spin network (F,7) we have 

(r,7)= + 



This implies in particular that the generating series Sr of Definition 14.11 satisfies Sr G Z[[ze,e G £^(F)]] 
proving one part of Theorem [S] 

Proof, (of Lemma [11. ip We will calculate (F,7)^ by writing out the contraction of tensors into the single 
contributions of basis elements in the various tensor products. Given a choice of permutations we place tensor 
products of basis elements e^^ (g) • • • (8) e^^ on every edge. Such a term contributes if some basis vectors on 
the same strand are unequal or the same basis element occurs twice in a tensor product corresponding to an 
edge. Otherwise the contribution is always (_i)#crossings^ 

To enumerate the non-zero states in an orderly fashion including the signs of their contributions we 
proceed in three stages. First we choose which sets of arcs get labeled by the same base vector. Then we 
assign base vectors to the loops in the curves. Finally we count the number of such labeled states. 

In the first stage we choose the collections of all loops that are labeled by the same base vector. For 
non-zero states these collections are necessarily 2-regular subgraphs of F, also known as curves. We will view 
curves both as sets of the loops they contain and as subgraphs of F. This first stage of the calculation is 
summarized in the concept of a curve configuration. To every curve c G Cr we assign the number of distinct 
base vectors it occurs with. The curve configuration C determines the total coloring of the spin network 7. 
However there may be many curve configurations yielding the same coloring 7. 

Note that all states with the same curve configuration contribute the same sign. This is because of 
planarity of the (multiple) annulus neighborhood of the curve: two loops in the same curve can only cross 
an even number of times. The common sign is determined by the cyclic ordering of the graph and equals 
(— 1)'('^^. This follows directly from a count of the number of intersections of such states. 

In the second stage we choose base vectors and we assign them to the loops such that all the loops with 
the same base vector constitute exactly one copy of a curve. First we need to choose which base vectors to 
use. There are (|^|) ways to do this. Then we need distribute these over the various curves, assigning £(c) 

base vectors to the curve c in no particular order. There are ('^') ways to do this. 

In the last stage we need to count how many arrangements of loops are possible given a curve configuration 
with base vectors at every loop. As was shown in |KLj there are T\ ways to do this. Their argument does 
not make use of planarity but for completeness sake we still reproduce the argument here. If we ignore the 
labelling by base vectors we will count exactly of the possible states. However if we also choose to ignore 
the fact that in contributing states the loops cannot self intersect, then counting possibly zero states will 
multiply our count by the same factor C\. To see why, note that to modify a zero state takes one permutation 
per curve. 

Recall that X! corresponds to a choice of a permutation on the arcs on each side of every vertex in the 
ribbon graph. We will show that there are exactly X! (possibly zero) states corresponding to a given curve 
configuration. First choose one such a state S. We will show that all others can be reached by inserting a 
unique permutation at every side of every vertex, with the appropriate modification. At any rate there cannot 
be more than X! such states, because arcs belonging to one side of a vertex can never be permuted to become 
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part of another side, because that would destroy the curve configuration. Even inserting a permutation 
at the side of a vertex may destroy the curve configuration, but in this case this can be compensated as 
follows. Choose an edge adjacent to the chosen side of the vertex and insert the inverse of the part of the 
permutation swapping arcs belonging to different curves. Exactly because of the fact that this compensating 
permutation involves only arcs corresponding to unequal curves this does not interfere with any other inserted 
permutations at the sides of vertices or their compensations. 
Summarizing, we get 

and rescaling we get 

■N\{\C\ 

Finally setting N = —2 provides a general formula for spin network evaluation as given in the lemma. □ 



To find a generating function for these evaluations, we first introduce a more general generating function. 
Definition 11.2. Define the curve generating Junction CGr of a graph F by 

c cec 

where CGr is a function of |Cr| variables, one for each curve in F. The exponents of these variables encode 
a curve configuration, and the evaluation (£) of a curve configuration C is given by 

(^>-(-ir<^'+|^|(i^i + i)(';^') 

The spin generating function Sr can be recovered from the curve generating function CGr by substituting 
Wc = IleGc ^e, becausc the color of an edge is the sum of the number of curves (with multiplicity) containing 
that edge. 

Proof. (Of Theorem [5]) Our goal now is to write CGr and hence Sr as a rational function. As a first step 
we note that if we drop the term (— l)'*^'''^ then we can sum the resulting series using the binomial series and 
the multinomial theorem to get 

E(-i)'^'(i^i + V' = (i+^. + .'..+^c.)- ^ ■ ■ ■ = 

Using this together with the Fourier expansion of (— f )'^('^' viewed as a function on the group {C2)'^^ as 

XcCr 

we find 

CGr = Y.{-iy^^hcw^^ E ^xY.i-^f'^bcw^ - E 

XcCr C XcCr 

where Qxiwi, . . . ,wcr) = Qi^i, ■ ■ ■ ,wcr) and Vi = ex{ct)wi. □ 

12. Challenges and future directions 

In this section we list some challenges and future directions. Recall from Definition 12.31 that every spin 
network (F,7) (and consequently every ribbon graph) has a spectral radius pr.7 and a number field Kr.-y- 
Our first problem implies the existence of an upper bound of the spectral radius in terms of F alone. This 
is phrased more conveniently in terms of the unitary evaluation of Definition [ 



Problem 2. Show that the unitary evaluation of a spin network (F, 7) satisfies 

l(r,7>''l<i 



26 



STAVROS GAROUFALIDIS AND ROLAND VAN DER VEEN 



This problem may be proven using unitarity and locality in a way similar to the proof that the Reshetikhin- 
Turaev invariants of a closed 3-manifold grow at most polynomially with respect to the level; see [Ga51 
Thm.2.2]. 

Recall the spectral radius of a cubic ribbon graph from Definition l2.3[ it is a nonnegative algebraic number. 
Our next problem is a version of the Volume Conjecture for classical spin networks with all edges colored by 
2. 

Problem 3. Show that the spectral radius of a cubic ribbon graph with 2n vertices equals to S'^". 

Theorem|3]confirms this problem for the following ribbon graphs: the Q, the tetrahedron, the 3-faccd drum 
and more generally for the infinite family of ribbon graphs of Remark 13.61 There is convincing numerical 
evidence for the above problem in the case of the i^a ^ and the Cube. 

Recall the generating series Fr^j{z) of Equation (jl]). Problem [3] motivates the next problem. 

Problem 4. Give a geometric meaning to the finitely many singularities of the generating series Fr^.y{z) for 
a spin network (r,7). 

At the moment, we do not know the answer to Problem |4] for the Cube with all edges colored by 2. 
Numerical evidence suggests that there are 12 equally spaced singularities in the boundary of the disk of 
convergence of -Fcube,2(^)- 

The next problem is a computational challenge to all the known asymptotic methods, and shows their 
practical limitations. 

Problem 5. Compute the asymptotics of the evaluation (A'3^3,2n) (given explicitly in Proposition [63| and 
(Cube,2n). 

Problem 6. Give a geometric interpretation of the number field Kr^-y of a spin network (r,7). 

In Theorem m we identified the Stokes constants of the 6j-symbol with geometric quantities of the asso- 
ciated metric tetrahedron. 

Problem 7. Give a geometric meaning to the Stokes constants of the sequence {r,n'y). 

The next problem is formulated by several examples involving the WZ method, discussed in Section 1101 
A positive answer can give an independent proof of Theorem [H 

Problem 8. Prove that for every coloring 7 of the tetrahedron spin network ( lt\,j), the sequence ( 1 , ^17) 
satisfies a second order recursion relation with coefficients polynomials in n. Can you compute the coefficients 
of this recursion from 7 alone? 

Our last problem concerns an integrality property of quantum spin networks. Our paper discusses four 
normalizations (r,7)^, (r,7), (r,7)^ and (r,7)^ for the evaluation of a classical spin network (r,7); see 
Definitions 11.11 fTT2l I6.2l and l9.4l With the notation of these definitions, these normalizations are related by 

(57) (r,7)^ = x!(r,7), (r,7)^ = ||(r,7), (r,7)^ = ^(r,7). 

On the other hand, if (r,7) is a quantum spin network, the Kauffman bracket evaluation (r,7)^ makes 
sense, and in general it is a rational function of the bracket parameter A; see Definition 16.21 and |KLj . The 
following problem is a first step towards the categorification of the evaluation of a quantum spin network. 

Problem 9. Show that for every quantum spin network (r,7), the evaluation {T £\/X\ lies in Z[A^^]. 

One can show that the above problem holds for the O and the tetrahedron spin networks, or more generally 
for the class of quantum spin networks of Remark 13.61 

Let us end this section with a remark. The main results of our paper can be extended to evaluations of 
spin networks corresponding to higher rank Lie groups. This will be discussed in a later publication. 
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